Terahertz magneto-spectroscopy of transient plasmas in semiconductors 
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Using synchronized near-infrared (NIR) and terahertz (THz) lasers, we have performed picosecond 
time-resolved THz spectroscopy of transient carriers in semiconductors. Specifically, we measured 
the temporal evolution of THz transmission and reflectivity after NIR excitation. We systematically 
investigated transient carrier relaxation in GaAs and InSb with varying NIR intensities and magnetic 
fields. Using this information, we were able to determine the evolution of the THz absorption to 
study the dynamics of photocreated carriers. We developed a theory based on a Drude conductivity 
with time- dependent density and density-dependent scattering lifetime, which successfully reproduced 
the observed plasma dynamics. Detailed comparison between experimental and theoretical results 
revealed a linear dependence of the scattering frequency on density, which suggests that electron- 
electron scattering is the dominant scattering mechanism for determining the scattering time. In 
InSb, plasma dynamics was dramatically modified by the application of a magnetic field, showing 
rich magneto-reflection spectra, while GaAs did not show any significant magnetic field dependence. 
We attribute this to the small effective masses of the carriers in InSb compared to GaAs, which 
made the plasma, cyclotron, and photon energies all comparable in the density, magnetic field, and 
wavelength ranges of the current study. 

PACS numbers: 78.20.-e, 78.20.Jq, 42.50.Md, 78.30.Fs, 78.47.-|-p 



I. INTRODUCTION 



The advent of long-wavelength coherent sources, such 
as free-electron lasers (FELs)J]l pacametric generators 
with difljerence frequency mixing,tl tecahertz (THz) 
antennas J3 and quantum cascade lasersci has brought 
an entirely new class of opportunities to study low- 
energy phenomena in solid state systems in the time 
domain and/or high- intensity regimes. In particular, 
far-infrared (FIR) / THz pulses can directly probe low- 
energy dynamics in bulk and quantum-cc)jD£n£d_aemicon- 
ductors, e.g., cyclotron resonauG|e--(CR),ffll3'ffiEji}itjj£^^aJ 
transitionSj-of shallow donoralilila and excijxtfis, 
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phonons,E£rEZl and intcrsubband transitions .E§ll 
In addition, small photon en^xeies enhance the pon- 
deromotive potential energyllSEj while precluding in- 
terband absorption and sample damage, leading to 
the possibility oLfxtreme nonlinear optical behavior in 
semiconductors E3'Ej 

In this paper we describe results of our study of the 
THz properties of photogenerated transient plasmas in 
semiconductors using a synchronized short-pulse THz — 
near-infrared (NIR) laser system with picosecond time 
resolution, both in the absence and presence of an ex- 
ternal magnetic field. By simultaneously monitoring the 
temporal evolution of the transmission and reflection of 
a THz probe pulse after NIR excitation, we carried out a 
dynamical study of the Drude conductivity of transient 
plasmas. More specifically, we were able to directly deter- 
mine the density and scattering lifetime of photocreated 
transient carriers as functions of time, i.e., n(t) and T{t). 

Another unique aspect of this technique lies in the fact 



that intraband FIR/THz spectroscopy is independent of 
whether the states involved are interband- active or not, 
thus providing a rare opportunity to directly probe non- 
radiative (or "dark") states. Dynamics involving such 
states are not observable with conventional interband 
transient spectroscopies, e.g., time-resolved photolumi- 
nescence spectroscopy. This unique ability makes it a 
powerful tool for providing insight into how optically cre- 
ated nonequilibrium electron-hole pairs lose their excess 
energies while relaxing toward the band edge through 
various scattering and thermalizatioft processes before 
eventually recombining to luminesce.E3 

In our previous work,§0 we described the first demon- 
stration of picosecond time-resolved cyclotron resonance 
(TRCR) of photogenerated transient carriers by moni- 
toring THz absorption as a function of magnetic field 
at fixed time delays between the NIR pump pulse and 
THz probe pulse. In the present paper, we systemati- 
cally investigated carrier relaxation in InSb and GaAs at 
different NIR intensities and magnetic fields. Our cal- 
culations based on a Drude conductivity with n(t) and 
T(i) successfully reproduced the main observed features. 
By fitting the theoretical reflectivity versus time to the 
experimental curves, we found a linear dependence of the 
scattering rate, 1/t, on the carrier depaity, n. This sug- 
gests that electron-electron scatteringEj is the main fac- 
tor in determining the carrier scattering lifetime inside 
the plasma in our density and time delay ranges. The 
application of a magnetic held in InSb resulted in dra- 
matic modiflcations in plasma dynamics whereas GaAs 
did not show any strong magnetic field dependence up 
to 8 Tesla. This behavior of InSb can be attributed to 
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FIG. 1: Schematic diagram of the experimental setup used for 
time-resolved two-color (NIR-THz) spectroscopy of transient 
plasmas in semiconductors. 



the small effective masses of its carriers, which made the 
plasma, cyclotron, and photon energies all comparable, 
under our experimental conditions, and their subtle in- 
terplay led to the observed rich plasma dynamics. Our 
theoretical simulations are in good qualitative agreement 
with the observations, supporting this explanation. 



II. EXPERIMENTAL METHODS 

The NIR laser source was a Ti:Sapphire laser seeding 
a regenerative amplifier. The amplifier produced intense 
NIR (Anir ~ 800 nm) pulses with pulse duration of ~ 200 
fs and pulse energies as high as ~ 1 mJ at variable repe- 
tition rates up to 1 kHz. The Stanford FELcIl produced 
pulses (in macropulses, as described below) with wave- 
lengths continuously tunable in the midinfrared (MIR) 
(3-15 nm) and FIR/THz (IjV^O /^m) with pulse dura- 
tions ranging from 0.6 to 2 psES and pulse energies as high 
as ~ 1 /iJ. In the present study, the FEL wavelength was 
mostly fixed to Athz = 42 ^m (or vthz = 7.1 x 10^^ Hz 
or huJTHz — 29.5 meV). 

A schematic diagram of our experimental setup for 
the two-color spectroscopy experiments is illustrated in 
Fig. ^ The NIR output of the Ti:Sapphire system was 
directed through a computer controlled variable delay 
stage, after which it was spatially overlapped with the 
THz beam from the FEL using a Pellicle plate. The two 
beams were thus made collinear as they were focused onto 
the sample using a parabolic mirror. The NIR pulse ex- 
cited nonequilibrium carriers across the band gap of the 
sample, which then absorbed a fraction of the incident 
THz probe pulse. The transmitted and reflected THz 



beams were then recollimated and directed to liquid ^He- 
cooled Ge:Ga photoconductive detectors. The THz out- 
put of the FEL was a pulse train of 10 Hz " macropulses" 5 
ms in duration. These macropulses each contained many 
(^ 60,000) ~ 1 ps duration " micropulses" separated by 
84.6 ns, corresponding to a repetition rate of 11.8 MHz. 
The micropulse-to-micropulse energy fluctuations were 
factored out using a THz reference detector (not shown 
in Fig. |l]) before the sample. The Ti:Sapphire oscilla- 
tor was locked to the seventh harmonic of this repetition 
rate, i.e., 82.6 MHz. Our synchronization electronics al- 
lowed us to select a single NIR pulse per FEL macropulse. 
A combination of the optical delay stage and electronic 
delays in the synchronization allowed for selective delays 
from to 84.6 ns with a picosecond resolution. With this 
arrangement, we were able to compare the intensities of 
the transmitted and reflected THz pulses before and af- 
ter the NIR pump pulse. The amounts of photoinduced 
change in THz transmission and reflection were recorded 
as functions of time delay. We then deflned photoinduced 
absorption as 
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where R and T are the reflectivity and transmissivity, 
respectively, which are functions of time delay, and Rq 
and To are their equilibrium values, i.e., before the arrival 
of the NIR pump pulse. For InSb, we determined, by 
Fourier transform infrared spectroscopy, that Rq = 18% 
(cf. theoretical value of - 20 %) and Tq = 9% for 42 ^m 
radiation. 

The InSb sample was undoped and had an electron 
density of 8.0 x 10^^ cm~'^ and a mobihty of 8.3 x 
10^ cm^V^^s^^ at 78 K. The GaAs sample was semi- 
insulating, with excess arsenic. We wedged both samples 
by ~ 3° to avoid multiple-reflection interference effects, 
and polished down to ~ 150 /Ltm, still much thicker than 
the absorption depths of both InSb and GaAs at 800 
nm. This has an impact on the carrier dynamics, as 
signiflcant carrier diffusion into the sample occurs (as 
discussed below). The sample was placed inside a 9 
T/1.5 K horizontal-bore split-coil magnet system (Ox- 
ford Instruments Spectromag 4000) with Sapphire cold 
windows and polypropylene room temperature windows. 
The sample was tilted 45° with respect to the magnetic 
field, B, which was parallel to both laser beams. 



III. EXPERIMENTAL RESULTS 
A. Power Dependence 

Typical zero-magnetic-ficld data for GaAs (left: a-c) 
and InSb (right: d-f) are shown in Fig. |^. The transmis- 
sion, refiection, and absorption of the THz probe beam 
are plotted against time delay. The wavelength of the 
THz probe was 42 ^m (z^thz = 7.1 x 10^^ Hz, Tiwthz = 
29.5 meV) and the sample temperature was 1.5 K. Each 
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FIG. 2: Temporal evolution of the THz transmission, reflec- 
tion and absorption for different NIR pump intensities in the 
GaAs (left: a-c) and InSb (right: d-f) samples. The wave- 
length of the THz probe was Athz = 42 fim (or vtuz = 7.1 x 
10^^ Ifz or fiLUTHz ~ 29.5 meV) and the sample temperature 
was 1.5 K. 
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FIG. 3: Photoinduced reflectivity change as a function of NIR 
pump intensity for GaAs for flxed time delays of ps, 418 ps, 
and 719 ps. Lines connecting data points are guides to an 
eye. It can be seen that with increasing NIR intensity the 
THz reflectivity can either increase or decrease, depending on 
the time delay. We can also see, for time delays of 418 ps 
and 719 ps, that the reflectivity initially decreases and then 
increases monotonically. 



panel shows multiple traces corresponding to different 
NIR intensities, with the maximum NIR fluence at the 
sample estimated to be ~ 4 mJ/cm^. In both the GaAs 
and InSb samples, the photogenerated carriers cause an 
abrupt drop (rise) in the THz transmission (reflection) at 
timing zero. For example, at the maximum NIR inten- 
sity, the transmission drop is 100% (complete trans- 
mission quenching) in GaAs and ~ 70% in InSb. The 
photoinduced absorption curves shown in Figs. ||(c) and 
y(f) were obtained from the measured transmission and 
reflection curves using Eqn. ([^). 

The subsequent recovery of these photoinduced abrupt 
changes depends critically on the relaxation properties of 
the sample under study. It is clear from Fig. g that there 
are significant differences between GaAs and InSb. In 
general, the GaAs sample shows smooth and monotonic 
temporal evolution throughout the entire time range pre- 
sented here (0—800 ps), whereas the InSb sample shows 
much more complicated behavior, exhibiting dynamic 
changes within the first 400 ps. The decay of the 
transmission (or absorption) change in InSb is far from 
monotonic, clearly showing multiple componets at high 
NIR intensities. Its reflection dynamics are even more in- 
triguing, exhibiting a sign change (positive to negative) 
at a certain time delay, which sensitively depends on the 
NIR pump intensity. As discussed in Section IV, we can 
explain these dramatic differences between the two sys- 
tems in terms of the importance of Auger processes, well 



known non-radiative carrier recombination especially im- 
portant in narrow gap semiconductors like InSbEj The 
much shorter absorption depth of InSb as compared to 
GaAs also favors the importance of Auger recombination. 
The sensitivity to the total carrier density, not the den- 
sity of the interband-active (or radiative) carriers alone, 
distinguishes the current spectroscopic technique from 
conventional transient spectroscopies based on interband 
transitions. 

The value of photoinduced reflectivity change is not 
a simple function of time delay or NIR pump intensity. 
Both its magnitude and sign depend on these parameters 
in a complicated manner. This is true even for GaAs, for 
which the reflection curves look smooth and monotonic 
in Fig. ||(b). To illustrate this point, we plotted the 
photoinduced reflectivity as a function of NIR intensity 
for three different time delays in Fig. |^. Here the data 
were taken by varying the intensity of the NIR pump 
pulse while the time delay was kept constant. We can see 
that with increasing NIR intensity the THz refiectivity 
can either increase or decrease, depending on the value 
of time delay. Also, if the time delay is large enough, 
we see that the reflectivity initially decreases and then 
increases with the NIR intensity. As we will see, the sign 
of the photoindiuced reflectivity change is governed by 
the interplay between the plasma frequency, ujp (cx ^/n) , 
and the THz photon frequency, wthz- 
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FIG. 4: Measured photo-induced reflection versus time delay 
at various fixed magnetic fields for (a) high (~4 mj/cm^) 
and (b) low (~15 /xJ/cm^) NIR excitations in InSb at 1.5 K. 
Traces are vertically offset for clarity. 



B. Magneto-Plasma Reflection 

Since GaAs did not exhibit any magnetic field depen- 
dence in our accessible magnetic field range, here we con- 
centrate on the data from InSb. Figures ^(a) and ^(b) 
show photoinduced reflectivity versus time delay for InSb 
at various fixed magnetic fields for high (~ 4 mJ/cm^) 
and low (~ 15 /iJ/cm^) excitations, respectively, at a 
temperature of 1.5 K. The data exhibit qualitatively dif- 
ferent behavior under these different excitation condi- 
tions. All the traces in Fig. ||(a) show significant pho- 
toinduced enhancement in reflectivity immediately after 
photoexcitation. This is because the plasma energy of the 
photocrcated carriers initially exceeds the photon energy 
of the THz probe (29.5 meV or 7.1 THz). The estimated 
initial carrier density in the high-excitation case is of the 
order of ~ 10^^ cm~^, which corresponds to a plasma 
energy of ~ 250 meV (or 60 THz). 

The low-excitation data at i? = T, on the other hand, 
shows a reflectivity drop due to the created carriers, as 
shown in Fig. ^(b). However, this completely opposite 
behavior very quickly disappears as we increase the mag- 
netic fleld from to 1.5 T. A small peak emerges in reflec- 
tivity, which grows in intensity with increasing magnetic 
field, reaches a maximum at ~ 3 T, stays roughly con- 
stant up to ~ 5 T, and finally goes away at higher mag- 
netic field. This very interesting behavior is explainable 
in terms o£-the tuning of the plasma energy by the mag- 
netic fieldo or the interplay between the plasma energy 



and the cyclotron energy, as discussed in more detail in 
Section IV.B. 



IV. DISCUSSION 
A. Plasma Dynamics 

We developed a theory for the reflectivity of a tran- 
sient plasma as a function of carrier density (n) , scatter- 
ing time (r), and magnetic field (B). After being created 
by the NIR laser beam, the carrier population decreases 
on a picosecond time scale. While there are a number of 
possible decay mechanisms, the dominant one, especially 
in narrow gap semiconductors at high densities (which 
is the case-Jos. in this study), is known as Auger 
relaxationm-BEl In this non-radiative decay mecha- 
nism, an electron and a hole recombine and the resulting 
energy is transferred to a third carrier. In this decay 
mechanism, the carrier density decreases in a character- 
istic way: 



dn 



(2) 



where C2 = 7.5 x_10~^ cm^/s is a reported Auger co- 
efficient for InSbEj In our analyses we used a modified 
Auger coefficient = KC2 in order to account for other 
decay processes, e.g., carrier diffusion and radiative re- 
combination, which are not explicitly taken into account 
in our model. The density evolution n(t) is then given 
by: 



n{t) = 



1 



C2t + l/n(0) 



t > 



(3) 



where n(0) denotes the initial density of photocrcated 
electron- hole pairs. 

The scattering rate is taken to be a power law 
function of the carrier density n, with coefficient a, ex- 
ponent f3, and a small fixed offset « 0.1 ps used to 
account for density- independent scattering mechanisms: 



(4) 



The dielectric function of a semiconductor after ultra- 
short pulse excitation of a nonequilibrium plasma is given 
by: 



1 -I- Cphiuj) 



Lu{uj — i/T[t)) 



(5) 



Here Coo = 15.68 is the dielectric constant for InSb, ujp 
= [4:TTn{t)e^ /meeoo]^^^ is the plasma frequency {nie — 
0.014mo). The phonon contribution Eph is calculated as: 



ep/i(w) 



(6) 
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FIG. 5: Gray scale image of the reflectivity of a plasma cre- 
ated in InSb as a function of carrier density n and scattering 
time T for a THz probe beam with a photon enenrgy of hui = 
29.5 meV. 
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FIG. 6: Experimental (a) and calculated (b) photoinduced 
reflectivity of InSb as a function of time delay for selected 
NIR pump intensities. Traces are vertically offset for clarity. 



where hcuL — 23.6 meV {Hujt — 22.2 meV) is the energy 
of the longitudinal (transverse) optical phonon and F = 
0.35 meV. For a fixed THz photon energy Tiwthz = 29.5 
meV (42 ^m, 7.1 THz), eph w - 0.17 - 0.005i is the 
material parameter. 

At 45 degrees incidence, the reflection coefficient 
can be calculated using a Fresnel formula: 



cos(7r/4) ~ J - sin2(7r/4) 

ru. (t) = \ . (7) 

e„ cos(7r/4) + -y/ - sin2(7r/4) 

In Fig. H, the reflectivity R — r^r^ for Tiwthz = 29.5 meV 
(t^THz = 7.1 THz) is plotted as a function of scattering 
lifetime, r, and carrier density, n, obtained through Eqns. 
(|), (|), and d), in the ranges of r = IQ-^^ - 10"" sec 
and n = 10^^ — 10^^ cm~^. The temporal evolution of 
the normalized photoinduced reflectivity due to transient 
plasmas is then computed as: 



m 

Ro 



= 1 



ARjt) _ rUtyjt) 
Ro r„(oo)r* (cx)) 



(8) 



which can be directly compared with experimental data. 
Here Rq — ri^{(X))r^{(X)) « 0.198 is the reflectivity in 
equilibrium. 

In Fig. |6| we present our experimental data (panel a) 
along with the fits (panel b) obtained using Eqn. (||). We 
used the initial rise of the refiectivity, i.e., AR{0)/Rq, to 
determine the initial carrier density, as well as the pa- 
rameters a and /? in Eqn. (^. The initial density of 
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FIG. 7: Temporal evolution of the density (a) and scattering 
rate (b) in InSb for selected NIR pump intensities. 



photocreated carriers n{0) scales linearly with the pump 
intensity, as expected. We also found that most fits result 
in a value of (3 very close to 1, [Suggesting the importance 
of electron-electron scatteringrJ The resulting expression 
for the scattering rate is as follows: 1/r « 2.1 x 10~^ n 
[cm^-^] + l/Ti [cf. Eqn. Using this expression for 

all NIR intensities, wc then performed the time-evolution 
fits assuming the Auger-like decay of the carrier popula- 
tion [cf. Eqn. (||)]. The multiplicative factor k showed a 
rather weak, but systematic, decrease from 17 to ~ 7 
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FIG. 8: Magneto-plasma reflectivity calculated by Eqn. (|n 
for a THz beam with huj = 29.5 meV as a function of plasma 
frequency (ujp) and cyclotron frequency (wc). One can see 
that the reflectivity can be larger or smaller than the equilib- 
rium value in a complicated way, depending on the values of 
carrier density and magnetic field. The following parameters 
are used: lot = 10.7, Hul = 23.6 meV, Hlot ~ 22.2 meV, and 
r = 0.35 meV. [cf. Eqns. (|) and (|).] 
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FIG. 9: Measured (a) and calculated (b) time-resolved photo- 
induced refiection in InSb at a range of fixed magnetic fields. 
Traces are vertically offset for clarity. 



with increasing NIR intensity. The excellent agreement 
seen in Fig. o makes us believe that the carrier diffu- 
sion must be completed in less than ^ 10 ps under our 
experimental conditions and its neglect in our model is 
justified. In Fig. we present the temporal evolution 
of the carrier density (a) and the scattering rate (b) cal- 
culated using Eqns. and (||) using the parameters 
derived from the fits. 



B. Magneto-plasma Dynamics 

When an external magnetic field is applied, Eqn. (|^) 
is modified, for CR active and inactive (— ) polariza- 
tions, respectively, into: 



1 + ephiio) - 



lu{lo ztujc — i/rit)) 



(9) 



where Uc — eB/meC is the cyclotron frequency. Corre- 
spondingly, Eqn. (0) changes to: 



cos(7r/4) - ^ e^^± - sin^(7r/4) 



e^,± cos(7r/4) -I- y'e^,± - sin^(7r/4) 



We then calculate the reflectivity, R, for a linearly- 
polarized THz beam with frequency u by taking an av- 
erage between the CR active and inactive circular polar- 
ization states: 



R 



This is plotted in Fig. g for a THz beam with Ucu = 
29.5 meV as a function of LOp and lOc- Here the value of 
ujT is fixed at 10.7, and Ulul = 23.6 meV, hcur = 22.2 
meV, and T = 0.35 meV are used for Eqn. Figure 
^ can be used for qualitatively predicting THz magneto- 
reflection dynamics by treating the Up axis as the time 
delay (since ujp oc y/n(t)) and the Uc axis as the applied 
magnetic field (since lUc (x B). Most importantly, one can 
clearly see from this figure that the instantaneous value 
of reflectivity of a transient plasma can become larger 
or smaller than the equilibrium value in a complicated 
manner, depending on the instantaneous values of the 
carrier density and magnetic field. 

The temporal evolution of the normalized photoin- 
duced reflectivity for a linearly-polarized THz probe with 
frequency lo is calculated as: 



m 



= 1 



Ai? ^ ^c.,+ (^)C+(t)+r^,-(OC^_(t) 
Rq r^,+ (oo)r* _^(c») -f ri^,_(cx))r* _(oo) ' 

(12) 

which can now be compared with experimental data. 

(10) Fi gure H illustrates (a) experimental and (b) theoretical 
magneto-plasma reflection dynamics. At zero magnetic 
fleld, the created plasma decreases the reflectivity. This 
behavior is due to the fact that the plasma energy of the 
created carriers is lower than the THz photon energy. 
When this is the case, an increased carrier density trans- 
lates into a decreased refractive index, and hence, a de- 
creased reflectivity. Even a small applied magnetic field 

(11) can dramatically change this behavior. A small peak ap- 
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pears at zero delay, whose magnetude rapidly increases 
with increasing magnetic field up to ~ 3 T, saturates, 
and gradually decreases. This intriguing behavior can be 
explained in terms of tuning of the plasma edge by B, 
i.e., the interplay between ujp and tOc, and is successfully 
reproduced by our theoretical calculation [Fig. ^(b)] . 

V. SUMMARY 

In summary, we performed picosecond two-color (NIR 
and THz) time-resolved spectroscopy on GaAs and InSb. 
We simultaneously monitored the dynamics of THz 
transmission and reflection while we varied the strength 
of the applied magnetic field and the time delay between 
the NIR and THz pulses. We found that the photoin- 
duced reflectivity dynamics are drastically affected by 
the magnetic fields. These results demonstrate the power 
and usefulness of this FIR technique for investigating the 



dynamics of nonequilibrium carriers in semiconductors at 
very low energy scales. 
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